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We formulate a holographic description of effects of disorder in conformal field 
^ ' theories based on the replica method and the AdS/CFT correspondence. Starting 

■ with n copies of conformal field theories, randomness with a gaussian distribution 

■ is described by a deformation of double trace operators. After computing physical 
^ ' quantities, we take the re — > limit at the final step. We compute correlation 

' functions in the disordered systems by using the holographic replica method as well 

as the formulation in the conformal field theory. We find examples where disorder 
changes drastically the scaling of two point functions. The renormalization group 
flow of the effective central charge in our disordered systems is also discussed. 
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1 Introduction 



Classical statistical or quantum many-body systems can be studied in the presence of 
spatial inhomogeneity. Indeed, it is quite rare that a perfectly clean system is realized in 
experiments, and hence having a good understanding of effects of randomness/impurities 
is to some extent necessary. Besides such a practical motivation, disorder by itself or com- 
bined effects of disorder and interactions can give rise to rich phenomena, which deserve 
studies in their own right. To name a few, complex behaviors in spin glass systems such 
as the Ising model with random bonds (random ferromagnetic interactions) or random 
magnetic field [1], or Anderson localization of electronic systems in the absence/presence 
of electron-electron interactions [21 [3] have been discussed. 

When the amount of disorder is small or disorder is (marginally) irrelevant in the 
renormalization group (RG) sense, effects of disorder can be studied perturbatively around 
a clean critical point. On the other hand, randomness is not necessarily small and it 
can drive the system to a new type of disorder-dominated critical point, called random 
critical point. Examples include a multicritical point in the random bond Ising model 
P,ll], the integer quantum Hall plateau transition [5], and a possibility of metal-insulator 
transitions in (2+l)-dimensional correlated electron systems [6]. These putative critical 
points are beyond perturbative or mean-field treatments and understanding the nature of 
these random critical points has remained to this date as a major challenge in condensed 
matter physics. 

It is the purpose of this paper to discuss quenched disordered systems in the framework 
of AdS/CFT correspondence j^. Effects of disorder have been studied as a perturbation to 
weakly coupled and exactly solvable (conformal) field theories, in terms of the perturbation 
theory or perturbative RG, assuming the disorder strength is small (see, e.g., [HI El [THl 
ITTl [12]). On the other hand, if we apply the AdS/CFT correspondence, it might be 
possible to solve strongly disordered problems because AdS/CFT is a strong/weak duality. 
In AdS/CFT setups, CFTs are typically non-abelian gauge theories [7| or the critical 
0{N) vector model [13]. Even though their precise relations to real condensed matter 
systems are not clear at present, recently there have been several hopes and circumstantial 
evidences that the AdS / CFT can capture essential features of condensed matter systems, 
such as the electrical and thermal transport [Hj, the quantum Hall plateau transition [T5] . 
the superconductivity [l6], the entanglement entropy [17], the scale invariant theories with 
non-trivial dynamical exponents [18], and so on. In the presence of weak randomness a 
holographic analysis for disordered systems was given in |19j . 

Assuming that disorder configurations such as the location of impurities are distributed 
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according to some underlying probability distribution, quantities of our interest (e.g., local 
correlation functions) also fluctuate from different disorder configurations. We are thus 
forced to deal with the probability distribution of observables or in particular the first (or 
first few) moment (s) of the observables. At this point, it is important to emphasize that 
the average over disorder configurations is taken after we take the statistical mechanical 
ensemble average over spins, or the path integral over (quantum) field configurations and 
so on. This is a major difficulty in disordered systems since for a given distribution 
of disorder we do not have translation invariance, although correlation functions after 
quenched disorder averaging may be translation invariant. A standard tool to analyze 
effects of disorder is the replica method (or replica trick) (see, e.g., [20 11211110] and section 
[2] in the present paper). 

In this paper, we consider a CFT deformed by a certain operator with its coupling 
changing from position to position. We will show that with the replica method, a gen- 
eralization of the double trace interaction [221 123] can nicely describe the randomness. 
Based on this idea, we will formulate the holographic replica method with some exam- 
ples in section El In particular, we calculate the two point functions and show that their 
scaling behaviors importantly change due to the randomness. In section IH we will give a 
complementary field theoretic analysis of the same system and confirm that it agrees with 
the holographic result. In appendix |Al we will present a generalized holographic replica 
model and realize the order /disorder phase transition. 

In contrast to clean systems, the structure of the RG flows connecting random critical 
points is less understood in disordered systems. This is so since critical theories describing 
random critical points are expected to be non-unitary. This in turn means we cannot use, 
in two dimensions, say, the c-theorem by Zamolodchikov, to know the direction of the 
RG flow. In particular, when the supersymmetric disorder averaging is applicable, the 
central charge is always vanishing c = because of the cancellation among matter fields 
and ghosts. Gurarie [2^ [T2] introduced, with the use of the supersymmetric disorder av- 
eraging, the effective central charge (or Gurarie's b) and proposed to use it to measure the 
degrees of freedom of a disordered two-dimensional CFT. Later, the equivalent quantity 
in the replica method (which we call Ce// in this paper) was introduced by Cardy [25] . 
(See appendix [HI) It is an open problem if the effective central charge Cg// shows the irre- 
versible relation Cgj?y < along the RG flow, as the central charge in two-dimensional 
unitary CFTs does. In appendix [Bl we take a first step of the holographic calculation of 
Ceff and show that the c-theorem like relation Cgj?j < holds in examples where we can 
calculate Cg// straightforwardly. 
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2 Disordered Systems and the Replica Method 



Let us start by reviewing how the repHca method can be apphed to quenched dis- 
ordered systems (see, for a review, e.g. [201 Ell [ID])0- Our arguments cover any quan- 
tum field theory (QFT) and we simply represent the action of a d-dimensional QFT as 
Sq[(p] = J d'^x Co{(p) with abstractly expressing all fields as (p. We pick up a certain opera- 
tor in this theory and denote it by 0{x). We perturb this theory by adding an interaction 
of the form 

S = So+ I (fx g{x)0{x). (2.1) 



This defines a classical disordered system in d spacial dimensions. In a (i + 1-dimensional 
quantum system, a spatially inhomogeneous perturbation by disorder is given by 
J dtd'^x g{x)0{x,t), where t represents the (imaginary) time direction, and the disorder 
configuration g{x) depends only on the spacial coordinates x. Below we shall proceed in 
the classical setup, though we can extend to the latter case straightforwardly. 

A quenched disordered system (or random system) is such a system where the coupling 
g{x) is depending randomly on the spatial coordinates x. We assume that the randomness 
is distributed with a gaussian profile, i.e., its distribution functional is given by 

P[(?(x)] oce-^^'''"^^")' (2.2) 



with / > 0. The free energy in the disordered system logZ can be found by simply taking 
the gaussian average of log Z over disorder, where represents disorder averaging with 
respect to the distribution P[g{x)]. 

In disordered systems, we are interested in correlation functions averaged over disor- 
dered configurations. An important point is that this random average is not equivalent 
to taking the random average for the action itself (12.11) . Instead the averaged correlation 
function is expressed as follows; 



{0{x,)0{x2) ■ ■ ■ 0{xk)) = (2.3) 

^/D(/pe-^«['^]-/^'^9(^)°(^)C(xi)C(x2) ■ ■ ■ C(xfc)^ 



Dg{x)P[g{x)] 



The most important technical problem is how to deal with the factor of the inverse 
partition function. One way to overcome this difficulty is to introduce ghosts and represent 



^ There are two types of disorder: annealed disorder and quenched disorder. The former is the case 
where impurities arc in thermal equilibrium with main system, while in the latter, it is not. In this paper 
we only discuss the latter case. Refer to [10] for more details. 
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the inverse partition function by that of ghosts. This is called the supersymmetric method 
and at a random fixed point it is typically described by a non-unitary CFT with the 
vanishing total central charge (see, e.g., [26l [12l [9] ) . This procedure is quite useful since 
we can utilize the knowledge of quantum field theory technique. However, it has a severe 
disadvantage that presently it can be used only for hmited situations, where the original 
theory described by S'o is a free field theory. 

Instead of proceeding along this direction, we would like to resort to another method 
called the replica method, which can be applied to any quantum field theories (refer to, 
for example, [2D])- The method may be summarized as follows. First we introduce n 
copies of the QFT and denote the z-th copy as QFTj. Next we prepare n copies of the 
field in QFTj. Then we consider the path integral in the product of n QFTs, i.e. 
QFTi ® QFT2 ® ■ ■ ■ ® QFT„, given by 

/n 
Dg{x)P[g{x)]\{[D^,] e-^"-i^°[^^l-/'^'^^(^)S"-^»(^)Oi(xi)Oi(x2) ■ ■ ■ Oi(x,) 

\[[D^.^ e-^-^°[-'l+^/'^^KEIUc..W) 0,{x,)0,{x,) ■ ■ ■ 0,{x,l (2.4) 
1=1 

where Oi{x) is the operator 0{x) in the Hilbert space of QFTi. In the replica method, 
we identify the average (12.31) with the n — limit of (12.41) . The parameter n is originally 
a positive integer number, but we assume the analytical continuity with respect to n 
and take the limit n — ^ finally. Indeed the inverse partition function required in (12. 3p 
formally appears in this limit. 

Let us define the scaling dimension of 0{x) by A, then the deformation induced by 
the randomness is relevant or marginal if 

2A < d. (2.5) 

This condition is clear from power counting in the replicated action fl2.4p and is called 
Harris criterion. 

We consider large gauge theories or a 0{N) vector model as examples of QFT since 
they appear in AdS/CFT correspondence. The operator O that couples to disorder g{x) 
is then given by a single trace operator, and with the replica trick the disorder effect is 
expressed by the double trace deformations as in fl2.4l) . which have been studied in the 
context of AdS/CFT [221 123], [27|. If we are interested in strongly coupled (or strongly 
disordered) regime of field theories, then the path integral in (12. 4p would be quite difficult 
to compute. Since the AdS/CFT correspondence maps strongly coupled CFTs to weakly 
coupled gravity theories, it is desirable to establish dual description to analyze strongly 



4 



coupled disordered systems. In this paper we mainly focus on the planar limit of large N 
gauge theory, which is dual to the classical limit of the dual gravity theory. 



3 Disordered Systems via AdS/CFT 

Now we move on to the main part of this paper: formulation of the AdS / CFT corre- 
spondence for disordered systems. We consider a standard setup of AdS/CFT, where a 
d dimensional Euclidean CFT is dual to a gravity theory on c? + 1 dimensional Euclidean 
anti-de Sitter (AdS) space described by the metric 

ds^ = + (3.1) 

The dual CFT is proposed to live on the boundary located at z = 0. 

A spin-less operator O, which is typically a single trace operator in a gauge theory, 
is supposed to be dual to a scalar fielcti with a mass m in the bulk AdS space. The 
relation between the conformal dimension A of (9 and the mass m is given by the formula 



A± = ^±^m2 + ^. (3.2) 

For our purpose, we need a relevant operator in the replicated theory, and hence we 
require 02.51) . From this condition we pick up the smaller dimension A_ in fl3.2|] and we 
call it simply A below. The behavior of near the boundary z = looks like 

^{z, x) ~ z'^-^iaix) + 0{z')) + z^ (^1^ + 0(z^)) . (3.3) 

In the standard interpretation of AdS/CFT [31], a{x) is regarded as a source to the 
dual operator O, while I3{x) is its expectation value {0{x)) = /5(x). We demand the 
normalizability for the mode with (j) z^, which leads to a constraint on the range of A. 
Combining with (12.51) we have to choos^ 

d — 2 ^ d 

— < A < -. (3.4) 

The lower bound is known to be dual to the unitarity bound in the dual CFT |31j . 



^We will stick to the scalar field example just for simplicity. We can generalize this to, say, vector and 
spinor fields straightforwardly. 

•^In order for A to satisfy this condition, the mass of </> should be in the range of —cP/A < m? < 
(1 — cf)/A. In AdS space tachyonic modes are allowed due to the curvature, and the lower bound is 
known as Breitenlohner-Freedman bound [30] . 



5 



For the scalar field we require the regularity at z = oo, then a{x) and j3{x) are 
related as follows 

In the momentum space representation it can be expressed as 

f]{k) = G{k)a{k), (3.6) 

where G{k) is given bjfl 

Since we chose A_ instead of A+ as the conformal dimension of dual operator O, the term 
proportional to a{x) is less singular than the one proportional to /3(x) in the boundary 
limit as opposed to the usual cases. Despite this fact we can still treat the former as a 
source because a{x) and (3{x) are related to each other via a canonical transformation as 
first claimed by Klebanov and Witten 



3.1 Double Trace Deformation in AdS/CFT 

According to [23j (see also [27j) a multi-trace deformation J df^x W{0) can be incorpo- 
rated into the formulation of AdS / CFT correspondence. Here we focus on the deformation 
of a CFT by a double trace operator of the form (A > 0) 

S^nt = \ j rf'x(0(x))^ (3.9) 

as a preparation for the later analysis of disordered systems. Assigning the boundary 
behavior (13.31) and demanding the regularity at 2; — > 0, the scalar field is uniquely param- 
eterized by a{x) or (3{x). As mentioned before these two parameters can be exchanged 
by a Legendre transform, and it is useful for our purpose to use /5(x) instead of a{x). In- 
serting the scalar field into the kinetic term and partially integrating over the coordinate 
z, we can obtain the action in terms of a{x). Then the Legendre transform leads to 

S^^\j^>^^. (3.10, 



*We employed the formula 



'^""^^^ ^ r(A) ^ ■ (^-^^ 
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Notice that the result is expressed as a field theory on the boundary of AdS space. 

Recall that the expectation value of 0{x) corresponds to the variable l3{x). Then the 
total action S in the presence of the double trace deformation (13. 9p and a source to (p is 
expressed as 



The equation of motion for (3 leads to 

1 + XG{k) 



^Pik) (-^ + A ) Pi-k) + Pi-k)Jik) 



(3.1i: 



m + J{k)=0. (3.12) 

In other words, the boundary condition for (p is now changed into a + Xf3 + J = 0. From 
fl3.12p we can express 

SM = -i/A.Jw(^^)^(-*), (3.13) 

and in the end we obtain the two point function 

taking the derivatives with respect to the source J{k). 

Due to the deformation of the (marginally) relevant operator 03.91] . the conformal 
symmetry should be broken in the CFT side. Nevertheless, the vacuum solution in the 
AdS side remains trivial 0(x, z) = even in the presence of non- vanishing A. Hence the 
background is still AdS^+i, which implies the conformal symmetry of dual CFT. This is 
so since our calculation neglects backreaction in the gravity theory, which should be taken 
into account from one-loop order [32]. Notice that the two point function fl3.14p breaks 
the conformal invariant In the UV limit A; — > oo, it is approximated by 

{0{k)0{-k)) ~ G{k) ~ (3.15) 

while in the IR limit we find 

{Oik)0{-k)) ~ i - ~ const. + Oik'-'"^). (3.16) 

This shows that under the RG fiow the operator with the conformal dimension A (= A_) 
fiows to the one with the conformal dimension d — A{= A+) [231 1221 [33]. In the IR 
limit conformal invariance is recovered and the operator dual to the scalar field has the 
conformal dimension A_|_. This does not cause any problems since the modes for both A_ 
and A+ are normalizable due to our restriction (13. 4p . 



^At the planar limit generic correlation functions do not receive any contributions from the deformation 
of double trace operator. This fact can be understood by explicitly writing down the double line Fcynman 
diagrams. The correlation functions involving O are exceptions. 
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3.2 Holographic Replica Method 




by employing the replica methoqj. If we consider a CFT with the random interaction 
(12. II) . it is described by a deformation of the product of n copies of the CFT, CFTi ® 
CFT2 ® ■ ■ ■ ® CFT„, in the replica method. Assuming the gaussian distribution of disorder 
(12.21) . the theory is deformed by double trace operators as (12. 4p 



where we denote Oi as the operator O with its conformal dimension A in CFTj. Com- 
pared with the standard double trace deformation (13.91) . the sign of the coupling (13.171) is 
opposite (since / in (12.21) must be positive), a common feature of the replicated theory in 
disordered systems. Thus, one may worry that this interaction may cause an instability 
in this theory. If we wish, we can add a conventional double trace deformation (13.91) to 
the original theory together with the randomness (12. ip . In the replica method, it means 
that we consider the product of n CFTs with the following generalized interaction 



Later, however, we will see that the limit A = of two point functions is well-defined, 
and the two point functions show more interesting scaling behaviors than the case with 
A > 0. 

In the dual AdS side, the spacetime is defined by n copies of an AdS space. They are 
disconnected in the bulk and attached to the same boundary aX z = Q. As seen above, 
the deformation by the double trace operator changes the boundary condition for the dual 
field, therefore after the deformation by (13.181) the fields in copies of an AdS space interact 
with each other through the boundary conditions. Interestingly, this setup with multiple 
AdS spaces has recently been discussed in [351 EB] from a different motivation. Originally 
we have n gravitons dual to n stress-energy tensors conserved independently. After we 
put the interaction (I3.18p . only a combination of stress-energy tensors is conserved and 
hence there exists only one dual massless graviton as argued in [551 EB]- In the AdS 
gravity the other n — \ gravitons become massive due to the one-loop contribution, while 
in the dual CFT, the conformal dimension of other combinations of stress-energy tensors 

slightly different replica method has been applied to AdS/CFT in [17] to compute the entanglement 
entropy [34] , 




(3.17) 




(3.18) 
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deviates from A = d to d + 6{n) due to the interaction fl3.18p [25| [35| [36] . The standard 
bulk to boundary relation relates 6{n) to the mass of n — 1 gravitons as M| = d6{n). 

In the following, we calculate two point functions at the tree level of AdS gravity in 
this setup. We denote as the field corresponding to Oi, and define and l3i{x) 

from the asymptotic behaviors as in fl3.3p . Following the previous analysis, we obtain the 
deformed action as 



S[p, J]„, = j d^k 



Er=iA(fc)A(-fc) 

2G{k) 



■ (3-19) 



A 



Here we included the source terms only for fii and without losing generality, by taking 
the symmetry into account. Notice that when A > 0, this system (13.191) is stable in the 
n — > limit ^. Taking the limit of A ^ in the end, we define the theory with A = 0. 
From the equations of motion 6S/6j3i = 0, we find 

^ -fG{k)Mk) - (1 + (A + (1 - n)f)G{k))Mk) 

{l + XG{k)){l + {X-nf)G{k)) ^ 
^ -fG{k).h{k) - (1 + (A + (1 - n)f)G{k))Uk) 

{l + \G{km + {X-nf)G{k)) ^ 
R- R fG{kf{J,{k) + Uk)) 

(i + AG(A:))(l + (A-n/)G(A;))- ^""-^^^ 

From these equations we can express the action S as 

Qrn -_1 f .'iJ Gik) (1 + (A + /(I - n))Gik)) {Mk)M-k) + Mk)M-k)) 
27^% il + XG{k)){l + {X~nf)G{k)) 

2fG{k)\h{k)U-k) 



f3 21 ^ 

'l + XG{k)){l + {X-nf)G{k)))' ^ ■ ' 

While n is a positive integer in our starting expression, in (13.211) we are free to regard n 
as a continuous valuable, which is a crucial assumption for the replica method. 
Two point functions can be computed from (13.211) as 

in(i.\n( Gik)il + iX + f{l-n))Gik)) 

mk)0,i-k)U = (i + AG(A;))(l + (A-n/)G(fc))' ^^'^^^ 

^^^^')^^^-'))^ - a^XGikmit-nf)Gik)y (^-2^) 



^Indeed, if we redefine $o = X)i Pi: Pi ^ Pi ^ '12 j Pj' then the interaction terms are given by 
i(A — nf )$Q + ^ ^ ■ $f. Therefore, any A > would be enough to stabihze the saddle points for small n. 
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The subscript n implies that the two point functions are evaluated in the replicated theory 
with fixed n. As discussed in section [2l we take the limit n — > and finally find in the 
disordered system as 

[l + {f + X)G{k))G{k) 



{0{k)0{-k)) = {0,{k)0,{-k)) = (i + XGikW ^^'^^^ 



We can also compute {0{k)){0{—k)), which involves, when replicated, two distinct repli- 
cas, 



fG{k) 



2 



{0{k)){0{-k)) = {0,{k)02{-k)) = (3.25) 



Unlike {0{k)0{—k)) , {0{k)){0{~k)) is made non-zero solely because of disorder, and 
hence our analysis indeed describes a disordered phase. Correlation functions of this type 
have been used as an order parameter in spin glass theories [1]. 

When A > 0, we can observe from the two point function (13.241) that the operator 
O with conformal dimension A always fiows into the one with A+ = d — A in the IR 
limit k ^ 0. Therefore, the IR limit of the disordered system looks similar to the theory 
deformed by a single double trace operator (13.91 ^. In the appendix |Al we suggest that 
the random spin system, such as the random bond Ising model [1], is analogous to this 
case with A = /, where the order/disorder phase transition occurs. In the appendix [B| 
we analyze the effective central charge Cg// [25| [8] of this disordered CFT and calculate 
the difference between the effective central charge at the UV (trivial) fixed point (c^^) 
and at the IR fixed point (Cg^j) . We find that the inequality c^^j- — c^jj < always holds 
for any values of A > 0. This may support an analogue of c-theorem in the disordered 
system, though there has been no proof of the c-theorem for Ce// from the CFT side. 

At the special point A = 0, we can obtain a markedly different result. In the IR limit 
k ^ 0, the operator O with dimension A becomes an operator with dimension 2A — d/2, 
as seen from the behavior 

{0{k)0{-k)) ~ (3.26) 

The constraint (13.41) leads to (i/2 — 2 < 2 A — d/2 < d/2, and hence the lower bound of 
2 A — d/2 violates the unitarity bound. However, this may be fine since the disordered 
system is not a closed unitary system. The analysis of the effective central charge Cg// at 
A = seems to need a special care as discussed in the appendix [B] and we leave it as a 
future problem. Observe that the IR limit is not fully gapped, but rather the theory flows 
into another critical fleld theory. As it becomes clear from the RG analysis in the next 



Fixed points for more generic cases with n = 2 are analyzed in |36| and their resuh is consistent with 
ours. 
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section, the random fixed point corresponds to the infinite randomness hmit (/ oo). 
This behavior could be comparable with random vector potential models of Anderson 
localization [37] or random quantum spin chains [38] . 

In principle, we can extend our computations to higher point functions by including 
interaction terms in (13.211) . In the case of a cubic coupling, we can add to fl3.19p a cubic 
term like 

Scubic = j dVk2 Cih, h) f3{h) f3{h) - h). (3.27) 

where C{ki, /C2) is related to the three point function in momentum space for / = via 
the standard rule [28l [29] . 

It is also possible to add more general multi-trace interactions f d'^xW{0{x)) in the 
boundary CFT. Correspondingly, in the holographic description, we need to add the 
potential term J d'^xW{P{x)) to fl3.19p . which leads to the generalized action 

2 



S = K I d^xd'y^-'^'^'^^^'^y^ 



d'^X 



■i Eft 



X 



i=l 



(3.28) 



where K is a numerical constant. 

Similarly, we can consider more general probability distributions for disorder other 
than the gaussian white-noise. After the theory is replicated, it leads to a deformation 



/ d'x,--- 1 d\g{x,)---g{xpf J] a,(xi)---a. 

p=l ^' ai,...,ap 



(3.29) 



where g{x\) ■ ■ ■ g{xp) is the p-th cumulant of the probability distribution P[g{x)]. 



3.3 Examples of Disorder Systems in AdS/CFT 

3.3.1 Random M = 4 Super Yang-Mills 

As the first example, we can consider the disordered system of A/" = 4 super Yang-Mills 
theory in four dimensions by a random deformation (12. ip . We can take (9 to be a 1/2 BPS 
operator with A = 2 made of two transverse scalars, i.e. Oat = Tr[$"$*] — |(5"''Tr[$'^$'^]. 
In this case, the condition (12.51) is saturated and therefore the RG fiow is logarithmic. 

3.3.2 Random 0{N) Magnet 

Klebanov and Polyakov [1^ conjectured that the massless fields in 74^5*4 with even 
spins describe the singlet sector of the three-dimensional critical 0{N) vector model in 
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the large limit. By using the vector field with N components, the action of the 0{N) 
vector model is 

1 /• _ r ^ ^ \ ^ ^ J 

(3.30) 



S = ^ I d-'x 



There are two critical points; One is at A = 0, i.e., the free field theory. The other is at 
the end point of the RG fiow induced by the second term in (13.301) . which is interpreted as 
a double trace deformation by setting O = (p ■ (p in section 13.11 Notice that the dimension 
of = {(p ■ (pY is 2 A = 2 and thus it is relevant. In the IR fixed point, the conformal 
dimension of O changes into d — A = 2 and then O"^ becomes irrelevant. 

Now, starting from the 0{N) vector model, we introduce the randomness via the 
interaction J d^xg{x) (j){x)-(j){x). In the replica method, this disordered system is described 
by the n — > limit of the following system 



S = ]^l d'x 



n f A 

,i=l *ii=l i=l 



(3.31^ 



Its holographic description is precisely given by the model (I3.19P analyzed in the previous 
subsection. We then conclude that whenever A 7^ 0, disorder is innocuous in the IR limit. 
In particular, disorder is an irrelevant perturbation at the non-trivial fixed point of (13.301) . 
which agrees with the well-known fact for the 3D 0{N) magnets with N > 2. (See, for 
example, [Ml SO].) 



3.4 Comments on Replica Symmetry Breaking 

The replica symmetry is the symmetry under a permutation of fields among different 
replicas, such as Pi [3j {i 7^ j) in our setup. The previous analysis (I3.19P clearly 
preserves the replica symmetry as all of j3i are vanishing. The breaking of the replica 
symmetry typically occurs when there are many vacua in the replica theory. If this 
happens, the analysis becomes more non-trivial because the definition of order parameter 
(or mean field) of randomness gets complicated. A famous such example is the problem 
of the spin-glasses (see e.g. P). 

Indeed, we can find examples where the replica symmetry is spontaneously broken in 
our holographic setup as follows. If we are interested in the IR limit, we can drop off the 
first term in (I3.19P and it is straightforward to study the vacuum structure of this system 
from the potential terms. By choosing W{13) appropriatel}|§, we can realize situations 



For example, we can realize this situation if we assume '^^i'^-* is a oscillating function so that its 



d/3 

average is increasing and '^^jf'^ = has only a single solution /? = 0. 
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where the potential W{l3i) has only the trivial vacu80, i.e., /3j = 0, while the potential 
with the term proportional to / has multiple vacua with /?j 7^ 



4 Field Theory Analysis in the Planar Limit 

In this section, we show that the two point functions fl3.24p and fl3.25p obtained from 
our holographic method can be reproduced from field theoretic calculations in the planar 
limit. This confirms the validity of our formulation of the holographic replica method. 
Before considering the replicated case, we reproduce from the field theory calculations the 
two point function (13.141) in the case of the deformation of fl3.9p 

S^nt = \j d^X^^U^), $pert(a;) = \{0{x)f . (4.1) 

In order to consider the renormalization of A and (9, we need to compute beta function j3x 
and anomalous dimension 70 . In the large limit, non-trivial contributions come only 
from the two point function {0{x)0{Q)) = with v = tt"*^/^ ^ ' those 

from higher point functions are suppressed (see, e.g., [23]). Therefore, divergent terms 
arise only from the following operator product expansions (OPEs) as 

2?? y 
$pert(a;)$pert(0) ~ ^^$pert(0), <l>pert (o^) 0(0) ~ Y^O{0). (4.2) 

From these OPEs we can obtain the beta functions following a standard analysis of quan- 
tum field theoriei"] 

-~Xik) = f3-, = {2A-dyX{k) + Cm?, io = A + h{k), (4.3) 



where we redefine A = Notice that even though we included only the 

leading order corrections to A, the result is exact in the large N limit. The beta function 
rSl) leads to 



Since the RG flow equation is given by 

^ {Oik)0{-k)) = 0, (4.5) 



d ^ d , ^ 



d\n\k\ dX 



^"^Here we include metastable vacua into our definition of vacua. 

^^For an extensive analysis of RG flows in the n = 2 case, refer to the third paper in 
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the two point function is 

(O(fc)O(-A:)) =Cexp p'\\r.\k'\{d - 2^o{k'))^ = YTJ^Y ^^'^^ 

The coefficient C is fixed such that {0{k)0{-k)) = G{k) for Aq = 0. 

Next let us introduce furthermore disorder (12. ip with the gaussian distribution (12.21) . 
In the rephca method the disordered system can be represented by introducing n copies of 
the CFT with Oi{x) (z = 1, 2, ■ ■ ■ , n) whose OPEs are Oi{x)Oj{0) ~ 6ijv/\x\'^^. Moreover 
we deform the n copies of the CFT by (13.181) 



S.. 



int 



^pert ("^y 



-/ j rf'^X$pert(x) + A j rf'^X^pert(x), (4.7) 
\i=l / 1=1 



As in the previous analysis, we now reproduce the two point functions fl3.24p and (13.250 . 
From the OPE coefficients, we read off the beta functions aJ^ 

j^^Rk) = = (2A - d)f{k) - n{f{k)f + 2f{k)~X{k), 
d 

d\n\k 
We find, therefore. 



.\(k) =p-^ = (2A - dfXik) + (A(A;))2. (4.9) 



- _ AW ^ (d-2A)(/,-V«) ^ ("-^AjA.^ 

n + Ao-n/o {kl'^-^^ + Xq ^ ' 

For computing the anomalous dimensions it is convenient to rotate the operators as 



— Y^Oiix), dj{x) = 0,{x) - - VOi(x), (4.11) 



n ^ — ' n 

i=l i=l 



with j = 1, ■ ■ ■ , ra, according to the irreducible representation of symmetric group 5„. 
Here the number of independent operators does not change since Yl]j=i ~ 0- ^^^^ 
normalization we have do{x)do{0) ~ v/\x\^^ and di{x)dj{0) ~ {6,j - i)tV|a;|2^. The 
anomalous dimensions of new operators are obtained just like before as 

1- n ~ 1~ 



When specialized to the case of the random 0{N) magnet discussed in section 13.3. 2[ the beta 
functions are consistent with the known results to the leading order in N. See, for example, [5^ HO] . 
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and solving the RG flow equation we have 



{Xo-nfo)G{ky \ ^ y/n nJl + \oG{ky 

(4.13) 

Rotating the operators back again, now we can reproduce the previous results of the two 
point functions (13.221) and (I3.23p . therefore after taking n — > limit we have (I3.24p and 

( Km . 

The method with the RG flow equation might be a standard way to compute cor- 
relation functions, but there is another way in a field theory viewpoint via a Hubbard- 
Stratonovich transformation first considered by Gubser and Klebanov [33] . The partition 
function we consider can be rewritten as 

Z][J]= (det-^y (^det^y J Dgix)p_Da,ix)x (4.14) 



where the subscript suggests that the correlation function is computed without pertur- 
bation, i.e. with / = A = 0. Integration over ai reproduces the double trace deformation 
f / d'^x J2iOf. Notice that the non-trivial contribution arises only from the two point 
function {Oi{x)Oj{0)) = SijG{x) (= 5jjf/|a;p^) in the large N limit. From this observa- 
tion we find 

Z][J] = (det-^y (det^) ' J Dg{x)f[Da^ix)x (4.15) 
where Gg{x) = J df"yG{x — y)g{y). Since this expression is gaussian with respect to g 



and cxj, we can easily integrate them out. The result io 

Z][J] = (1 + AG)-'^(1 + (A - n/)G)-5x (4.16) 

with 

Q = — Q' = — (4.17) 
1 + AG l-nfQ 

Taking derivatives with respect to Jj twice, we obtain the two point functions fl3.22p and 
fl3.23p and again reproduce fl3.24p and (I3.25P in the n ^ limit. 



^■^The prefactor corresponds to the 0{N^) corrections to the partition function, and it can be used to 
compute the shift of central charge along the RG flow. See, for example, [33] and appendix IB] 
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5 Conclusions and Discussions 



In this paper, we have studied the quenched disordered systems with arbitrary strength 
of disorder by applying the AdS/CFT correspondence. We formulate the holographic 
replica method by employing the setup of the double trace deformation in AdS/CFT. We 
have calculated the two point functions in our disordered system in the planar limit from 
both the AdS and CFT sides and got the same result. We found that the scaling of the 
two point functions evolves non-trivially under the RG flow. Especially, if we fine tune 
the parameter to A = 0, then the two point functions in the IR limit show remarkably 
new behavior. In the generic case A > 0, the IR limit is essentially the same as the fixed 
point obtained by the standard double trace deformation. As in appendix [B] it is possible 
to analyze the effective central charge Ce// using the AdS/CFT correspondence. There 
we observe that it decreases under the RG flow between two fixed points when A > 0. 
A more thorough analysis of this would be very intriguing as no proof of the c-theorem 
about Ceff has been known so far. 

It is important to remember that the subtle limit n ^ of the replica method does not 
cause any problems in our examples. We also pointed out that the replica symmetry may 
be broken if we consider particular deformations of the CFT. It may also be intriguing to 
note that the limit n — offers us a formal way to construct AdS duals to non-unitary 
CFTs with the central charge c = 0. 

Even though our holographic formalism of the replica method covers quantum dis- 
ordered systems, our explicit computations are performed only for classical disordered 
systems. Thus here we would like to mention the application of our holographic replica 
method to quantum disordered systems. Since the random coupling g{x) does not depend 
on the time t, the frequency u of the field f3i in (13.191) should be vanishing. Thus the re- 
sults fl3.24p and (13.251) remain the same only when u = 0; Otherwise we will get the same 
result as the one in the pure system, i.e., with / = 0. This triviality is because we are 
taking the planar limit in the presence of the multi-trace interactions. In order to obtain 
a;-dependent results, we need to go beyond the tree level analysis in the gravity theory. 
A similar situation occurs when we are considering two point functions of other opera- 
tors, such as the currents J'^{t,x). At the tree level, we do not find any /-dependence in 
both classical and quantum disorder systems. These problems of one loop analysis clearly 
deserve future studies. 

In addition to the replica method, we know another method which also enables us to 
study disordered systems, called the supersymmetric method, as mentioned before. In 
order to deal with the randomness in a theory with scalar fields and Dirac fermions we 
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add their superpartners (or ghosts). In this procedure, the bosonic global symmetry, such 
as 0{N), becomes its supergroup extension, such as OSp{N\N). It may be interesting to 
find a similar method for Yang-Mills gauge theories. Naively, one may think that we can 
replace the bosonic gauge group U{N) with the supergroup U{N\M) (see [H] for gauge 
theories with supergroup gauge symmetries). However, we can easily see that this leads 
to what we precisely want only at the vanishing coupling qym = 0. Therefore we need a 
modification or another idea. 
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A Order/Disorder Phase Transition 

In condensed matter systems, randomness often competes with ordering tendencies. 
For example, for a spin system (e.g. the Ising model) with critical temperature Tc, the 
system is in an ordered phase (e.g. ferromagnetic phase) for T < Tc. Even if we introduce 
a small amount of disorder the system will still be in the ordered phase. However, if the 
randomness becomes strong enough, a phase transition will occur and the system will 
evolve into a disordered phase. This is a standard story, for example, in the ferromagnets 
in the presence of random magnetic field (see e.g. jlOj). 

Such competition can be described in our holographic approach. In the spin system 
examples, we regard O as the spin operator a. Then an important point is that we omit 
the i = j terms in the interaction terms — ^ T^,- ^ OjOj (see fl3.18p ) since the operator 
: : does not exist in spin systems like the Ising modeo Another point is that when 
T < Tc and / = 0, the system is in a ordered phase. As an example we express this 
with the spontaneous breaking of the Z2 symmetry (3 — s> —/3 by adding the standard wine 

-'^^For example, if we consider the Ising model, the OPE of the spin operators is <t ■ a = I + £. This 
OPE produces the energy density operator £, but this just shifts the critical temperature. In contrast, 
if we consider the large N 0{N) vector model, then the operator : : exists and we have to keep it as 
in (|3.19p . Even in the latter case, we can add the interaction term as in p.lSp with A = / and realize a 
system without the i = j terms. 
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bottle potential W{/3) = + X(3^. Notice that > when T < Tc, while = 

at T = Tc- In summary, we reach the following action of (3 



S[(3]^ = njd'^k 



In the n — >■ limit, the term due to the randomness behaves like and it competes with 
the spontaneous breaking of the Z2 symmetry. Therefore we can realize the disorder /order 
phase transition when we increase the randomness parameter / in our holographic de- 
scription. 



(A.l) 



Assuming that the vacuum preserves the replica symmetry, i.e., /^i = = ••• = /?„ = 
we can rewrite (lA.ip as follows 



(A.2) 



B Effective Central Charges in Disordered Systems 

We can measure the degrees of freedom of a given CFT in even dimensions by calcu- 
lating the central charge. Let us suppose that the replica theory has a non-trivial fixed 
point. Then we can define the central charge of this new CFT as 

c(n) = nco + Ac(n), (B.l) 

for any fixed n. Here cq is the central charge of the original CFT without any deformations. 
For example, in the A/" = 4 U{N) gauge theory case it is given by cq = The term 

Ac{n) is due to the deformations and comes from one-loop order corrections 0{N^). In 
the random system, we can still define so called the effective central charge Ce// (see 

dc{n 



Ceff = 



(B.2) 

n=0 



dn 

Physically, this central charge is equal to the coefficient of the OPE of energy stress tensors 
at the random fixed point [25] 

) OC Ceff. (B.3) 

In a usual unitary CFT, there is a c-theorem, which states that the central charge is 
decreasing under the RG fiow. However, in disordered systems, no c-theorem is known 
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for Ceff even if the replicated theory with fixed n is unitary, because finally we need to 
take the formal n ^ procedure. 

In our setup of AdS/CFT, the order 0{N^) correction to the central charge is pro- 
portional to the one from the one-loop effective potential \4_ioop = — |TrAds(~n + "^^) 
of the scalar field. Since our formulation of holographic replica model is closely related 
to the case with the double trace deformation (13.91) . let us first review the case analyzed 
in [32]. In the case with the double trace deformation (13. 9p . the A-dependence of Vdoubie 
comes from the boundary condition for the field (j){x, z) imposed at the boundary z = 0. 
For a{x) and (3{x) in (13. 3p . we assign 

a{x) = -X(3{x), (B.4) 

which can be obtained from the relations (13.121) and (13. 6p . The one-loop vacuum energy 
was explicitly calculated in [32] as 

^douMe(A) = - 2.-.,./2r(rf/2)i^^+i X ^' r(z>)r(i-z>) I (^-^^ 

where i/ = | — A and A = 2^'" p||^^| Xz'^'^. Here z is the radial coordinate and R is the 
radius of AdSd+i- In the above expression, we subtracted the divergent piece which does 
not depend on /, and which would actually be canceled by other contributions due to 
the super symmetry. From this result we can show that Vdoubic(A = oo) < 0, and this is 
consistent with the c-theorem as it is proportional to cjr — cuv- 

Now we compute the one-loop vacuum energy for the product of n copies of the CFT 
with the deformation (13.181) . In this case the boundary conditions are modified as 

n 
i=l 

for j = 1, 2, ■ ■ -, n, or performing the redefinitions as in footnote [3, 

"0 = -{X - nf)po, oij = -Xpj. (B.7) 

Here we should notice that J2j = J2j f^j = by construction. Therefore, by comparing 
with the previous result, we find 

l^replica(/, A, u) = V^ioublc (A - uf) + {u - 1 ) ^double ( A ) . (B.8) 

Notice that the identity Vrepiica(/, A, n = 0) = holds as expected. What we are interested 
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in is the shift of the effective central charge Cg// (IB.2P under the RG flo\\o 

dAc{n 



Ac, 



eff 



dn 



= A(Klouble(A)-/-Kl'ouble(A)), (B.9) 
n=0 



where A is a positive constantl^i. 

We would like to calculate the difference between Cgf/ at the UV fixed point A = / = 
and at the IR fixed point. If we assume A > 0, then Ce// of the IR fixed point can be 
obtained by setting \ = oo (recall the analysis of the RG flow in section Hj). Then we 
immediately find 

Ce// - = A ■ Kiouble(oo) < 0. (B.IO) 

This seems to support the c-theorem like property for the effective central charge of our 
disordered CFT. However, we would like to stress again that the c-theorem for Cg// has 
not been proven at present and a counter example may be found in more generic cases. 
This value of Cgj?j — c^j^ (IB. 101) is actually the same as that of the standard double trace 
deformation [321 [33] even in the presence of the random perturbation. 

Moreover, if we restrict ourselves to the case with < / < A, then we can show 
that Ace// is always non-positive for any values of / and A from the expression ( IB. 51) . 
Therefore we conclude that Acg// is a monotonically decreasing function of / and A in 
this case. 

In the remaining special case A = 0, we would get naively the opposite result c^jj — 
^eff ^ However, this analysis is not reliable because Vdoubic {—nf) in ( IB. 81) is divergent 
for any values of n > 0. This subtlety arises since at A = the system of ( I3.19P becomes 
unstable as mentioned before. We may instead deform the theory by, for example, a 
quartic term so that it become stabilized. We would like to leave the analysis of c^jf 
at A = as a future problem. 
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